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Abstract 

We develop a Luttinger liquid theory of the Coulomb drag of persistent 
currents flowing in concentric mesoscopic rings, by incorporating non-linear 
corrections to the electron dispersion relation. We demonstrate that at low 
temperatures, interactions between electrons in different rings generate an 
additional phase and thus alter the period of Aharonov-Bohm oscillations. 
The resulting nondissipative drag depends strongly on the relative parity of 
the electron numbers. We also show that interactions set a new temperature 
scale below which the linear response theory does not apply at certain values 
of external flux. 
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During the last decade, persistent currents (PC's) in mesoscopic rings have attracted 
significant interest both theoretically [0-0 and experimentally 0-^. Much of this atten- 
tion was due to a large discrepancy between experimentally observed current amplitudes in 
disordered metallic rings 0J^, and theoretical predictions based on a single-particle picture 
0. Yet unresolved, this puzzle has generated a number of theoretical works on the role 
of electron-electron interactions in multi-channel disordered rings. 

At the same time, in clean single-channel rings the theory predicts that PC's at low 
temperatures should not be affected by interactions |1T0|,^,|6[| , and exhibit Aharonov-Bohm 
oscillations as a function of flux with the same period and amplitude as for non-interacting 
electrons. Results of a recent experiment on a single semiconductor ring with low number 
of channels are in agreement with these predictions . 

On the other hand, interactions should become essential in a system consisting of a pair 
of clean ID rings with different radii, placed concentrically, as is shown in Fig. 1. If the 
rings were isolated, PC's in each ring would oscillate with a period determined by its radius. 
Inter-ring interactions will change the oscillation pattern by causing Coulomb drag of PC's, 
as we show below. 

There are, in general, two physical mechanisms for the current drag. The flrst mecha- 
TT[] , which originates from "friction" between two subsystems caused by scattering of 



nism 



carriers in one subsystem by density fluctuations in the other, has been widely studied dur- 
ing recent years [T^, following experimental observation of the Coulomb drag [0. At low 
temperatures, the resulting transresistance behaves as (T^ In T in a disordered system) 
and vanishes then with decreasing T as the phase space for scattering shrinks. The second. 



nondissipative mechanism was pointed out by Rojo and Mahan |]T4|, and is based on the 
observation that the van der Waals interaction between two current carrying subsystems is 
modifled, resulting in a flnite current drag at zero temperature. 

The latter mechanism is, in fact, responsible for the PC drag in mesoscopic rings. Inter- 
ring interactions change the ground state energy of the system. As a result, PC in each 
ring, being a thermodynamic quantity, acquires a dependence on the flux penetrating the 
other ring. With decreasing temperature, this dependence should become sharper because 
of a singular ("sawtooth") shape of the PC as a function of flux at zero temperature. In the 
vicinity of a "tooth" , even a small perturbation may affect strongly the amplitude of the PC. 
For this reason, an adequate theory of the PC drag must be valid beyond linear response. 

In this paper we develop such a theory using the Luttinger liquid (LL) approach to PC 
in ID rings. This approach is based on Haldane's concept of topological excitations in flnite- 
size ID systems [^], and was recently extended by Loss to account for parity effects |10 



in the presence of external flux. However, in its standard formulation, the LL approach does 
not allow for current drag, which arises from the asymmetry between electrons and holes. 
With electron dispersion linearized, the asymmetry is lost and the electron and hole drags 
completely compensate each other. In LL formalism, the electron-hole symmetry manifests 
itself in a complete separation between zero modes, which carry PC, and bosonic flelds, which 
are responsible for inter-ring interactions. We derive an interaction-induced correction to the 
zero mode spectrum by incorporating the lowest order non-linear correction to the electron 
dispersion relation. Further, we show that at low temperatures, the entire effect of inter-ring 
interactions is to generate a phase which, in turn, leads to PC drag by changing the flux 
seen by the electrons. 
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We begin by recalling the LL approach to interacting (spinless) electrons in a ID ring with 
circumference Lj = 1, 2) threaded by magnetic flux 0j |1T5|.|^. By linearizing the spectrum 
near the Fermi points, the electrons are decomposed into right and left moving fermions with 
ipiaix) {a = ±) satisfying twisted boundary conditions ipia{x + Li) = (—1)^^ e'^'^'^^^'^° ^jJia{x) , 
where the sign factor reflects the dependence of the ground state on the parity of the particle 
number Ni and 0o is the flux quantum. The fermion operator can be expressed through 
bosonic fields as ipia{x) = (2eLj)"^/^e~*'^^°(''\ with ipia{x) = [6i{x) — aipi{x)]/2. The fields 
iPi{x) and 9i{x) can be presented as fl^ 



iPi{x) = (pij + Mi2-Kx/Li + (Pi{x), 

Oi{x) = + {Ji - 2$,)27r(x + L,/2)/L, + 



(1) 
(2) 



where Mi and Ji are fermion number and (topological) current operators, and 6iM and 
ipij are their conjugates: [9iM,Mi'] = [yjjj, Jj/] = 2i6ii'. For further convenience, we have 
incorporated both the flux and parity dependence into a single quantity $j, which is defined 
as $j = 0j/0o for odd, and = 4>i/4>o — 1/2 for Ni even. Owing to the boundary 
condition, the eigenvalues of Ji and Mj satisfy the selection rule that their sum is even. The 
periodic fields and 6i{x) have the usual representation in terms of bosonic creation 

and annihilation operators 



e^{x) 
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1/2 



e'^'^iaL+aikJ, 



Liki 



1/2 



sgn(A;,)e*"(al^^ -a,fcj. 



(3) 
(4) 



where fcj = 2np/Li {p being an integer), and the regularization factor exp{—e\ki\Li/27i) is 
implicit. In the absence of backscattering, ipi is related to the charge density (relative to the 
background density Ni/Li) by pi = 8x^1/211. The eigenvalues of Mj represent the numbers 
of extra electrons added to the ring while that of Ji are the excess of left over right moving 
fermions |[T5[] . 

In terms of the bosonic fields y^j, the (normal ordered) Hamiltonian of free fermions with 
linear dispersion reads 



(5) 



where Vi is the Fermi velocity and Ilj = Pi — ^i/ Li = OxOi/Air, Pi being the canonical 
momentum. 

In calculating PC's it is convenient to go to the Lagrangian formulation and present the 
partition function Z in terms of a functional integral. The free (Euclidian) action has the 
form 



'S'O = E / o [(^rV^i)^ + Vii^xV 

i JO Jo OTTVi 



i(^i + Si)dr'^i 



(6) 



where (3 is the inverse temperature and ip{T,x) has a decomposition (up to a constant) 
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iPi{T,x) = 27miT/P + 2nmix/Li + (/}j(r, x). (7) 

Here zero modes nj and rrii are winding numbers in r and x directions, respectively, and 6i in 
(D enforces the selection rule by taking values Si = (1/2) for rrii even (odd) Note that 
rrii coincides with the eigenvalue of Mj. In the following we assume the number of electrons 
in the rings fixed and restrict ourselves to the rrii = (with 5j = 0) sector. 

The partition function is given hj Z = J D{pe~^ , where S = Sq + Sint- Here the term 



1 ^ rP fLi rLj 

Sint = -^22 J dr J dx J dx'd^(piVij{x,x')d^>ipj, (8) 



with Vij being the Coulomb potential, describes intra- and inter-ring interactions, and the 
measure Dip includes the sum over zero modes together with the integral over periodic fields 
Tpj^. The PC is found by separating out the contribution from the zero-modes. With ipi{T, x) 
of the form (|^) (with rrii = 0), the latter are completely decoupled and the action takes the 
form 



'5 = Eh;F-27rm,$, +5, (9) 




where S = So[Tp] + Sint[^] is flux- independent and Tj = 2vi/'KLi is the temperature scale set 
by the finite size of each ring. Thus, the PC, J7i = f3^^d,j)- In Z, is carried by zero modes whose 
spectrum is unaltered by inter-ring interactions. This is a consequence of the linearization of 
the electron dispersion and the resulting quadratic form of the bosonic action 5*0. Physically, 
linearization lifts the electron-hole asymmetry, which is responsible for the current drag, so 
that PC in each ring is sustained by its own flux 0j. 

Let us now consider the lowest non-linear correction to the fermion dispersion relation. 



Following [|T5|, we take the correction to the linearized fermion Hamiltonian as 

2^ E £^ dx : 4M^li^^c.{x) : , (10) 



where rrie is electron mass. Transforming (|T0|) with the help of bosonic representation of 
ifjiai^x), we obtain 

Hni = E r • + 3^-¥'^(47rn,)' : , (11) 

which means that the corresponding bosonic theory is no longer free. Regarding Hni as a 
small term, we can obtain a correction to the action by adding (|ll|) to (^ and repeating the 
steps leading to the functional-integral representation of Z To first order in 1/m, the 
correction takes the form 



Sni = — J2 [ dr [ dx 

ASirme i Jo Jo 
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{d^^if ^d^(pi{driPif 

vf 



(12) 

Using (|^) (with mj = 0), we separate out the zero-mode contribution in Sni by writing 
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E 



iPnievf Jo 



dr / dx (f^drdx^Pi + S 



nil 



(13) 



where Sni depends only on fields Tp^, which are now coupled to the zero modes via the first 
term. Adding (|13]) to (P), we then perform the functional integral e~^^ = J Dlpe^^^^"' . 
Note that after the zero modes are separated out, the small term Sni can be neglected and 
the remaining gaussian integral iexplicitly evaluated. The resulting F{ni,n2) admits an 
expansion in terms of l/mg with the odd orders vanishing due to translational invariance. 
Thus, to the first non-vanishing order 



F(ni,n2) =Fo + ^F, 



ijTliTlj ^ 



(14) 



where Fq is and 0j independent. The second term in (0), which is to be combined 
with the first term in (P), represents the correction sought to the zero-mode spectrum. The 
diagonal coefficients Fa can be absorbed into Tj via renormalization of the Fermi velocities 
and do not play any role in the following Factorizing out the zero- mode part of the 



(15) 



partition function, Zq = e = e^^°Z, we finally arrive at 

n\n2 \ P ij i J 

where an = and the parameter au = P^Fi2, given by 



-1 



dZi / dZ2dr^dx^Di2{Zi,Z2)dr2dx2D2l{z2,Zi) 



(16) 



describes the coupling between zero modes in different rings. Here Zi = {Ti, Xi) and Di2{zi, Z2) 
is the nondiagonal part of the boson Green function, Dij{zi, Z2) = (^4(2:1)^^(^2)), calculated 
from the action S. Due to the azimuthal symmetry of the system, a convenient basis for 
Dij is given by angular-momentum eigenfunctions, in which the Fourier-transform of the 
inverse Green function reads 



4:7lViLi 



2 , 2/27rp 



+ 



LiLj 



(17) 



where Vij{p) is the Fourier-transform of the Coulomb interaction. For d Li the interaction 
is given by 



27r(p+ l/2)d 



(18) 



where k is the dielectric constant, Kq is the modified Bessel function, and we used the 
notation dij = d ioi i ^ j, and da = w, w being the width of the rings. The expression for 
ai2 then takes the form 



-TT 



ai2 



Y.u^y[D^2{u,p)]'. 



(19) 
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The frequency sum in (|T9D is straightforward and yields 



ai2 



dv 



12 



1 



(20) 



where 



t± = tanh 



Qiip) = jAip)±Bip), 



A{p) = [r(l + vu) + r-\l + V22)]/2, 



2e2 



B{p) = ^[r{l + Vu)-r~^{l + V22)]y^ + v 
"27r(p+ l/2)rf, " 



2 

12) 



(21) 

(22) 
(23) 



and r = (Ti/T2)^/^ characterizes the asymmetry between rings. Let us obtain an estimate 
of ai2 for identical rings, with Vi = v, Ti = Tq, and Lj = L. The analysis of (pOf) shows that 



for T <^ Til = Tq{L / d) ^ 1 + a\n.{2d / w) , where a = 2e^ / ktiv is dimensionless interaction 
constant, the coupling ai2 is temperature independent. Since for ^ L one has Td ^ Tq, 
and To is the crossover temperature above which the PC is exponentially suppressed, ai2 
can be regarded as a constant in the entire range of relevant temperatures. For d ^ w the 
screening is strong and ai2 is determined by the first order in the inter-ring interaction. The 
rhs of (|20) then gives 



a 



ai2 



32Tod^kl[l + a\n{2d/w)f/^'' 



(24) 



where fc^;' if the Fermi momentum. A d~^ dependence (without logarithmic factor) was 
obtained previously for infinite wires [|14|. 

Turning to the partition function (|15|), we observe that since the zero modes in different 
rings are now coupled, the PC in, say, ring 1, J7i = —dFo/dcpi, depends also on the flux 
through ring 2. Since the coupling ai2 is small, one could try to obtain a correction to 
the current by expanding the free energy to first order in ai2. It is important to realize, 
however, that at low temperatures such an expansion does not exist for all values of $j. In 
order to make this point clear, let us rewrite (pIS]) in a different form using Poisson's formula 
(omitting constant prefactor) 



Zo=J2 exp -Pj2(Pi ~ *i)cij(Pi - 



(25) 



PlP2 



where c = vr^a"^. For ai2 ^ an one has for diagonal elements Cn = n'^Ti which, in the 
absence of interactions, is the level spacing at the Fermi level. With such form of Zq, the 
PC is given by 



Jl($l,<l>2) = — 



Cii((pi) - $l) +C12((P2) - '^'2) 



(26) 



where (pi) stands for the average of pi calculated from the partition function (p5|). Note that 
$i in (^) takes values in the interval < < 1, and J'i is periodically continued outside 
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of this interval. For further analysis, it is convenient to obtain an equation for (pi). This 
can be done by making use of the following identity 



1 1 °° / 
{P^) = 2^2^ (tanh/3 

n=— oo \ 



^ Cij{pj - %) + Cii{n + 1/2 - Pi) 



(27) 



which can be readily derived by substituting Jacoby's product formula for the 6'3-function 
|T9| in place of the sum over variable pi in the double sum (^) (one notices, for example, 
that for C12 = the rhs of (|27| ) is p-independent and {pi) — 1 reduces to the log-derivative of 
the ^3-function). For low temperatures, T ca, the fluctuations of pi are suppressed and 
all moments factorize, (p"^) = (pi)"^, which allows one to replace pi by (pi) in the rhs of (pT]). 
It is also easy to see that for such temperatures, all terms in the sum with n 7^ cancel each 
other out (up to exponentially small corrections). Then (|27|) simplifies to the system 



iPi) 

{P2) 



/o[cn(l-2$i) + 2ci2((p2)-$2) 

/o[c22(l-2$2)+2Ci2((pi)-$l) 



(28a) 
(28b) 



where fo{x) = (e^^ + 1)~^ is the Fermi function. This system, together with (|2^) determines 
the PC at low temperatures, T -C Cjj = vr^Tj. In the absence of inter-ring coupling (ci2 = 0) 
one recovers the PC for an isolated ring. 



j;0($0 = 2/. /o 



Cud - 2$,; 



(29) 



where li = cu/ipQ = eVi/Li is the current amplitude. The current J'^ is a linear function 
of flux except in the interval $j — 1/2 ~ T/cu <^ 1 in which it rapidly changes from — /j 
to /j. Note that $j depends on the parity of the total number of electrons, iVj, resulting in 
diamagnetic (paramagnetic) current for Ni odd (even) [|10|,3- 

Turning the coupling on, PC's can be found by iterating the system (|28|). First, in the 
argument of the Fermi function we substitute (pi) expressed via J7i by neglecting the second 
term in ( pB]) (and in a similar expression for ^2)- Substituting (pSf ) back into (p6D, we finally 
obtain 



J,i^,,^2)=Ji°i^i-S,), 5i=6:72($i,$2), (30) 

with b = 0oCi2/2ciiC22 = — 0o«i2/27t^. Eq. (|30|), which is our main result, describes the mu- 
tual dependence (Coulomb drag) of PC's in coupled rings. It shows that electrons encircling 
the ring i acquire additional (Berry's) phase proportional to PC in the ring j. Thus, inter- 
ring interactions change the period of the Aharonov-Bohm oscillations. In particular, peak 
positions of PC get shifted by an amount Remarkably, the resulting oscillation pattern 
depends strongly on the relative parity of electron numbers Ni, as it can be seen in Fig. 2. 
We emphasize that in the vicinity of the peak, $j ~ 1/2, PC ( |5DD cannot be expanded in 
6i if the temperature is low enough: T ^ T* = C12. Similarly, J7j ($i, $2) in the expression 
for 6i cannot be replaced by J7j°($j) for $j 1/2. If $j is not too close to 1/2, one can 
neglect 6i in the argument of the Fermi function and the correction to PC takes the form 
SJ^i/Ii = '2b J j (note that b is positive). 
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Eq. (|30|) was obtained for temperatures lower than Tj. With increasing T, the PC 
amphtude gets damped due to the loss of phase coherence. Expanding the exponent in (^Sf) 
in terms of C12, the correction to the free energy can be written as 

(5Fo = 0ofeJi°(*i)J2°(*2), (31) 

with = (4/j/7r)(T/Ti)e"^/^' sin(27r$j) [g]. Correspondingly, the correction to PC for 

T ^ Tj is suppressed by an additional damping factor. 

Finally, let us address the experimental implications of our result. PC's with amplitude 
of about 4 nA were observed in a micron-size GaAlAs/GaAs ring with width w ~ 160 nm, 
and Fermi wavelength of electrons Af — 40 nm, which corresponds to a ~ 0.5 and Tq ~ 200 
mK. Assuming that two concentric rings were fabricated on the same plane, the estimate of 
ai2 for the above parameters and d ^ w yields a small phase shift 6 ~ 10~^. However, this 
small value is due to the rather large distance between rings in the plane. The situation will 
be significantly improved if two concentric rings were fabricated on parallel 2D layers. Since 
the interlayer distance can be as small as 3 nm this would allow one to decrease the 
actual separation between current carrying channels by 1-2 orders of magnitude. A reliable 
estimate of 012 for d -C w would require careful treatment of the screened interaction between 
the two closely spaced rings and depend on details of the structure. The ci-dependence ( pi]) 
suggests, however, that in this case one might expect a phase shift of several percent, which 
would be within resolution of experimental observations. 
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FIGURES 



FIG. 1. Schematic picture of two concentrically placed rings. 

FIG. 2. Example of PC in ring 1 at (j)2/4'i = 1-1 and Tj/T = 5.0 for different inter-ring 
couplings: cu/ca = (solid line), and cu/ca =0.1 with the same (dashed line) and opposite 
(dotted line) parities of electron numbers. 
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